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We investigate the time evolution of the quantum meson modes in the late time of 
chiral phase transition. In particular, it is shown that there exists a possible solution 
to the equation of motion for the quantum meson modes, which reveals a parametric 
resonance and/or resonance through forced oscillation induced by the small oscillation of 
the chiral condensate. After that, we demonstrate the unstable regions for the quantum 
meson modes in both the cases of a uniform and spatially expanding system. 



1. Introduction 

In the experience of ultra-relativistic high energy nucleus-nucleus collision, it is 
concerned that chiral symmetry is restored. However, when it is cooled down im- 
mediately, spontaneous symmetry breaking occurs again. When chiral condensate 
returns to actual vacuum, the amplitude of quantum pion modes with low mo- 
menta even in the late time of chiral phase transition may grow. This phenomena 
may connect to the parametric resonance. 

In this paper, we are going to introduce the other possibility which may exist 
to this mechanism, especially we would like to focus on the possibility by forced 
oscillation. To investigate the dynamical chiral phase transition, we use the vari- 
ational approa ch t o the 0(4) linear sigma model in the Gaussian wave functional 
approximation.f22l if^ jg shown that the amplification occurs in the mechanism of 
the resonance by forced oscillation as well as the parametric resonance induced by 
the small oscillation of the chiral condensate.^ 



2. TDVA with Gaussian wave functional 

In this section, we find the equations of motion for the chiral condensate and quan- 
tum fiuctuations. In order to treat the chiral condensate and the fluctuation modes 
around it self-consistently, we use the time-dependent variational approach (TDVA) 
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with a Gaussian wave functional based on a functional Schrodinger picture. We start 
with the following Hamiltonian, 

where a runs from to 3. The index indicates the sigma field and 1—3 indicate 

the pion fields. The small h means explicit chiral symmetry breaking term. There 
are three model parameters used here, i.e., to, A and h. The model parameters m^, 
A and h are taken so as to reproduce the pion mass = 138 MeV, the sigma 
meson mass Mq = 500 MeV and the pion decay constant /7r(= (fio) = 93 MeV in 
vacuum, expHcitly, = -(519.70 McV)2, A = 84.96 and h = (78.13 McV)^. 

We adopt the Gaussian wave functional as the trial wave functional in the 
framework of the functional Schrodinger picture. The time-dependence of varia- 
tional functions is governed by the time-dependent variational principle: 

5 j dt(i^^-H^=Q . (2) 

As the result, we obtain the following equations of motion: 

+ Ml{t) - -^^Mtf^ Mt) = h, (3) 

{d^ + k' + M!it)}u^it)=0. (4) 

These form a set of basic equations of motion for the chiral condensate and quantum 
meson fields. The ipo can be regarded as chiral condensate and the ui"' can be 
regarded as mode functions for quantum meson modes. 

We now demonstrate qualitatively the time evolution of the mean field (fa and 
of the quantum meson mode functions ui"^ in the case of a uniform system without 
spatially expansion. We assume that (po ^ and <fi = 0, with i = 1 — 3, that is, 
the chiral condensate points in the sigma direction. In the numerical calculation, 
we adopt the box normalization with length L in each direction. We then impose 
periodic boundary conditions for the fluctuation modes, namely, the allowed values 
of momenta are kx = {2tt / L)nx and so on, whore Ux is an integer. The fluctuation 
modes labeled by (n^, ny,nz) are included in each direction up to n"^ = n1 + riy + 
rj2 = 8^. This corresponds to the three momentum cutoff A ^ IGeV (more precisely, 
990 MeV), since we have adopted the collisional region as = (10 fm)^. Further, 
we assume that the fluctuation modes of the pion fields are all identical, and we 
denote them Ui = U2 = = u^- 

In Fig. 1, the time evolution of the chiral condensate is depicted. In order to avoid 
the complexity about problems of the initial conditions in relativistic heavy ion 
collisions, we only demonstrate the time evolution qualitatively with (^o{t = 0) = 40 
MeV and (fio{t = 0) = 0. The reason we adopted the above initial conditions is that 
we are interested in the behavior of the time evolution of the chiral condensate and 
quantum meson modes in the late time of the dynamical chiral phase transition 
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Fig. 1. The time evolution of the chiral condensate. The horizontal and vertical axes represent 
time and the value of the chiral condensation, respectively. 
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Fig. 2. The time evolution of the fluctuation mode with n = in the 7r-direction. 

in this paper. It is shown that the chiral condensate approaches the vacuum value 
with oscillation. 

In Fig. 2, it is shown that the amplitude of the quantum pion mode with n — 
is amplified. Similarly, the amplitude of the pion mode with n = 1 is also amplified 
weakly. However, amplification is not realized for modes with n > 2. By contrast, 
the sigma meson modes with n > are not amplified. 

Thus, we conclude that, in this parameterization, the amplitudes of the lowest 
and first excited pion modes are amplified, but other quantum meson modes are 
not. 

The time evolution of the chiral condensate with the expanding geometry is 
also calculated. It is seen that the damped oscillational behavior is realized more 
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quickly than in the case without spatial expansion. 

It should be noted here that the amplification occurs in the late time of the 
chiral phase transition in both the cases of no spatial expansion and one-dimensional 
spatial expansion. It is thus seen that, even when the chiral condensate oscillates 
around its vacuum value with small amplitude, there are amplified solutions of pion 
modes. The mechanism of the amplification is clarified in the next section. 

3. Late time of chiral phase transition 

In this section, we investigate the mechanism of the amplification for meson modes. 
It is shown that one of possible mechanisms is a forced oscillation, as well as para- 
metric resonance. 

Let us investigate the time-dependent solutions around the static configurations. 
The condensate and quantum meson modes can be expanded around the static 
solutions. 

<fo{t) ^(pQ + Sip{t) , 

"aW=^fW+<5«aW- (5) 

Here, we consider the late time of the chiral phase transition. Then, |(5<y5(<)| is small 
compared with the vacuum value, ipQ. Further, we assume that \Su^\ < With 
the above approximation, wc have following equation, 

Sulit') = F„cos(7i')e-'"°*' , (6) 

where a = a or n. Here, we introduce the dimensionless variables (cJq, , /iq, , Fq, ) . If 
Fa (-Pjr) is negHgible, then Eq. ([6]) is reduced to Mathieu's equation. In this case, the 
existence of the unstable solution for 6u^{t') may be expected. This phenomenon 
is well known as parametric resonance in classical mechanics. On the other hand, 
forced oscillation may be realized through the effect of (-FV), even if the model 
parameters do not allow for unstable regions for the parametric resonance. These 
phenomena are seen in the lowest and first excited pion modes, as demonstrated in 
the numerical calculations presented above. 

Let us search the the unstable regions for Eq. ([5]). We perform the time in- 
tegration over a certain interval for the magnitude of the quantum meson mode 
function, which means a kind of time average. Then, we compare the time integral 
at a certain time step with a proceeding one. If the value of this integral is larger 
than the last one for all time steps under consideration, then we decide that this 
mode is unstable. 

Figure. 3 shows the unstable region for quantum meson modes in the case of 
Fq. = 3.0. The vertical axis represents and the horizontal axis represents uj'^ha/2. 
The time integration is taken in 10 fm/c, in order to judge the stability of solution 
up to 100 fm/c. It is shown that the unstable regions become wider due to the 
effect of the forced oscillation and/or the beat induced by the oscillation of the 
chiral condensate, as compared with the parametric resonance only. 



— 2 +w^[l - /i„ cos(7t')] 
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Fig. 3. The unstable regions for the solutions in Eq. ^ with Fa = 3.0. 
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Fig. 4. The unstable regions for the solutions in Eq. (O with Fa = 3.0. 

In the case with spatially expanding geometry, the time interval for the time 
integration to determine the unstable regions in our approach is taken as the rather 
short time 4 fm/c, because the unstable regions disappear when there is a strong 
damping effect for a long time. We calculate the time evolution up to 40 fm/c. The 
damping effect induced by the spatial expansion is introduced in the equations of 
motion for the quantum meson modes as follows: 

i^-f^+,ul[l-h^coa{jt')e-^'']^du^^^{t') = F^cos{^t')e~^^^ . (7) 

In the case of Fa = 0, the result for the unstable regions is well known, because 
this case corresponds to the case of the Mathieu equation with a friction term. 
Here, in Fig. 4, we demonstrate the case of i^a =3 in Eq. (7). In this case, the 
forced oscillation and beat may realize. However, these effects only give a small 
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modification. The reason for this is that the friction has a strong effect in this 
model. As a result, the amplification of the quantum meson modes with the spatially 
expanding geometry is not strong, although the amplification appears. These results 
imply that the effect of the forced oscillation or the beat ceases to be effective for 
the spatially expanding geometry. 

4. Summary 

We have demonstrated that the amplitude of quantum pion modes are amplified 
even in the late time of the chiral phase transition through the mechanism of forced 
oscillation or the beat, as well as the parametric resonance, in the framework of the 
0(4) linear sigma model. Also, we have determined the parameter regions in which 
the unstable solutions for the quantum meson mode functions exist. Further, we 
have found that the effect of forced oscillation ceases due to the strong friction in 
the case with the spatially expanding geometry. 

Of course, it is necessary to know the initial conditions in order to judge whether 
the amplification of the quantum meson mode occurs or not, that is, whether the 
parameters are in the unstable region or not, in realistic relativistic heavy ion 
collisions. 
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